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ABSTRACT 

On July 4, 1997, after traveling close to 500 million 
km, the Pathfinder spacecraft successfully completed 
entry, descent, and landing, coming to rest on the sur- 
face of Mars just 27 km from its target point. In the present 
paper, the atmospheric entry and approach navigation 
activities required in support of this mission are discussed. 
In particular, the flight software parameter update and 
landing site prediction analyses performed by the Path- 
finder operations navigation team are described. A suite 
of simulation tools developed during Pathfinder’s design 
cycle, but extendible to Pathfinder operations, are also 
presented. Data regarding the accuracy of the primary 
parachute deployment algorithm is extracted from the 
Pathfinder flight data, demonstrating that this algorithm 
performed as predicted. The increased probability of 
mission success through the software parameter update 
process is discussed. This paper also demonstrates the 
importance of modeling atmospheric flight uncertainties 
in the estimation of an accurate landing site. With these 
atmospheric effects included, the final landed ellipse pre- 
diction differs from the post-flight determined landing 
site by less then 0.5 km in downtrack. 

NOMENCLATURE 

g deceleration constant, 1 Earth g = 9.806 m/s 2 
gi first deceleration measurement, nominally 5 g, g 
g 2 second deceleration measurement, nominally 
sampled 12 sec after gx, g 
t g0 computed time from g 2 point to parachute 
deployment, sec 

dtdg inverse of the nominal deceleration slope at 
gi, sec/g 

dgdt nominal deceleration slope at g 2 , g/sec 
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g 2 a minimum acceptable g 2 deceleration 
measurement, g 

g 2 b maximum acceptable g 2 deceleration 
measurement, g 

tgo2a maximum acceptable t go computation, sec 
t g0 2b minimum acceptable t g0 computation, sec 
y atmospheric entry flight-path angle, deg 

INTRODUCTION 

On July 4 1997, the Pathfinder spacecraft ushered 
in a new era of planetary exploration by successfully land- 
ing on the surface of Mars. The primary objective of the 
Pathfinder mission was to develop and demonstrate a low- 
cost, reliable system for landing on the surface of Mars. 1 
In addition to this engineering objective, a focused set of 
science investigations were performed and several tech- 
nology elements required for further exploration of Mars 
were demonstrated. 2 " 4 

In the present paper, the atmospheric entry and ap- 
proach navigation activities required to support 
Pathfinder’s entry, descent, and landing operations are 
discussed. These trajectory simulation analyses were 
applied by members of the Pathfinder operations navi- 
gation team to ensure successful parachute deployment 
and estimate the Pathfinder landing site. While the July 
4 entry sequence was autonomous, ground controllers 
had numerous opportunities to update onboard software 
parameters and adjust the atmospheric entry conditions 
during interplanetary cruise. 

The analysis tools used in entry operations were de- 
veloped in the spacecraft design and development phase 
and used to prescribe spacecraft test conditions prior to 
launch. The validity of these analyses is currently being 
analyzed through post-flight trajectory reconstruction. 5 
After a brief mission summary, the present paper con- 
tains an overview of the atmospheric analysis tools, a 
discussion of the parachute deployment algorithm devel- 
oped with these simulation analyses, and a description 



of their use by the Pathfinder operations navigation team 
in the days preceding the spacecraft’s July 4 entry, de- 
scent, and landing. Insight into the performance of the 
Pathfinder parachute deployment system based on a pre- 
liminary assessment of the flight data is also presented. 

MISSION OBJECTIVES 

The primary engineering objective of Pathfinder, 
acquisition and return of data on the entry, descent, and 
landing system as well as lander performance, was com- 
pleted a few hours after landing. Other objectives, 
including deployment of the Sojourner rover, were com- 
pleted on July 5, 1997. Figure 1 shows the lander on the 
surface of Mars on July 4, 1997; whereas, Fig. 2 is an 
image taken by the rover of the lander after its deploy- 
ment. While the primary mission duration for the rover 
and lander was 1 week and 1 month respectively, there 
is nothing to preclude longer operations. In fact, both 
systems have already surpassed this expectation, yield- 
ing a rich scientific harvest. 

The Pathfinder science payload is comprised of the 
Sojourner rover, an alpha-proton X-ray spectrometer 
mounted on Sojourner, a mast-mounted 360-degree ste- 
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Fig. 3: Mars Pathfinder science instruments . 

reo imager, and an atmospheric science/meteorology sta- 
tion. These instruments, depicted in Fig. 3, have 
permitted investigation of the landing site geology and 
surface morphology, the magnetic and mechanical prop- 
erties of the surface materials, determination of the 
atmospheric structure, and insight into the rotational and 
orbital dynamics of the planet (inferred from high-gain 
antenna tracking). 6 



Fig. 1 : Pathfinder spacecraft on the Mars surface, 
July 4, 1997. 



As shown in Fig. 4, the Pathfinder spacecraft was 
targeted for a landing within a 100 x 200 km ellipse in 
the Ares Vallis floodplain region of Mars centered at 
19.24 N latitude, 33.1 W longitude. This landing site is 
approximately 850 km southeast of the Viking 1 Lander. 
An Ares Vallis landing was selected for several reasons, 
including spacecraft and rover design constraints, entry, 
descent, and landing concerns, scientific potential, and 
safety. 7 Engineering factors which affected this choice 
of landing site included the desire for a low surface el- 





evation (to allow sufficient parachute deceleration time), 
high Earth and Sun elevation angles (for improved com- 
munication and power), and a relatively flat region free 
of surface hazards such as large craters, chasms, or knobs 
(for improved airbag performance). From a science per- 
spective, this site was deemed favorable since a variety 
of rock and soil types thought to be deposited by an an- 
cient catastrophic flood were expected. 

THE ROUTE TO MARS 

The Pathfinder spacecraft began its 7-month jour- 
ney to Mars with the launch of its Delta-II 7925 rocket 
from Cape Canaveral, Florida on December 4, 1996. To 
satisfy numerous mission constraints a Type I interplan- 
etary trajectory was selected with a heliocentric transfer 
angle of 155 deg. 8 Among these mission constraints, the 
requirement of a pre-dawn landing in the Northern hemi- 
sphere dictated a retrograde atmospheric entry at Mars 
(with the drawback of increasing the vehicle’s velocity 
relative to the atmosphere by approximately 0.43 km/s 
relative to a similar posigrade entry). During its interplan- 
etary cruise, a series of four trajectory correction maneu- 
vers were executed as shown in Fig. 5. 9 As discussed in 
a subsequent section, a contingency fifth trajectory cor- 
rection maneuver could have been performed at either 
entry -10 or entry -5 hours; however, when these deci- 
sion points approached (entry -1 1 and entry -7.5 hours), 
the flight team deemed this maneuver unnecessary. 



Fig. 5: Mars Pathfinder interplanetary trajectory . 

As shown in Fig. 5, the total transit time for the Path- 
finder spacecraft was 2 1 0 days in which the vehicle trav- 
eled approximately 5.0e+08 km along an ellipsoidal arc 
about the Sun. From an orbit determination analysis based 
on all available Deep Space Network tracking data at the 
time of entry, the best-estimated Mars atmospheric entry 
conditions were derived and are presented in Table 1 . In 
this table, both the inertial and relative velocity estimates 


are given. At entry, the spacecraft had a mass of 585.3 kg 
and a ballistic coefficient of 62.4 kg/m 2 . 

Table I: Mars Pathfinder best-estimated atmospheric 
entry conditions (10:00 am PDT, 7/4/97) 


Radius, km 3522.200 

Declination, deg 22.630 

West Longitude, deg 21.831 

Inertial velocity, km/s 7.264 

Inertial flight-path angle, deg -14.060 

Inertial azimuth angle, deg 253.148 

Relative velocity, km/s 7.479 

Relative flight-path angle, deg -13.649 

Relative azimuth angle, deg 253.675 


ENTRY, DESCENT, AND LANDING 

The Pathfinder spacecraft utilized a low-cost but 
complex entry, descent, and landing (EDL) strategy to 
survive its flight through the Mars atmosphere. As shown 
in Fig. 6, four deceleration mechanisms (aeroshell, para- 
chute, solid-rockets, and airbags) were used to slow the 
spacecraft from its interplanetary approach velocity (7.48 
km/s relative to the atmosphere) to its final velocity of 
zero. As shown in Fig. 7, entry, descent, and landing 
was initiated thirty minutes prior to encountering the 
atmospheric interface (defined at a radius of 3522.2 km) 
when the cruise-stage was jettisoned. The aeroshell en- 
countered a peak heat rate of approximately 1 00 W/cm 2 



Cruise-Stage Separation Entry and Descent 


Landing Surface Operations 

Fig. 6: Pathfinder entry, descent, and landing system. 
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M Cruise-stage separation (8500 km, 6100 m/s) 

™ Landing - 35 min 

Entry ( 1 33 km, 7479 m/s) 

^ Landing - 5 min 

Parachute deployment (9 km, 375 m/s) 
Landing - 171 sec 


Heat shield separation 
Landing - 149 sec 

f Lander separation/bridle deployment 
Landing - 1 30 sec 


Radar ground acquisition (1.5 km, 60-75 m/s) 
Landing - 25 sec 

Airbag inflation (300 m, 50-60 m/s) 

Landing - 10 sec 

Rocket ignition (50-70 m, 50-60 m/s) 
Landing - 6 sec 

Bridle cut (0-30 m, 0-25 m/s) 
Landing - 2 sec 



Deflation/petal latch tiring 
Landing + 15 min 


Final retraction 
Landing + 180 min 


Airbag retraction 
Landing + 115 min 


Fig. 7: Entry , descent , and landing sequence of events ( reconstructed from flight data). 


and a peak deceleration of 16 g approximately 70 and 
78 sec after encountering the Mars atmosphere. The para- 
chute was unfurled 171 sec past the entry interface. Re- 
lease of the forebody heatshield and extension of the 
bridle followed 22 and 41 sec after parachute deploy- 
ment. Once the spacecraft traveled down the bridle, the 
radar altimeter initiated the search for the Mars surface 
(at a surface altitude of approximately 1.5 km). This al- 
timeter data was then used to infer the descent rate and 
determine the appropriate time of airbag inflation (296 
sec past the atmospheric interface) and solid-rocket ig- 
nition (300 sec past the atmospheric interface). After the 
three solid-rocket motors were fired to delete the remain- 
ing vertical velocity, the bridle was cut (at an altitude of 
21 m above the surface). The spacecraft then fell to the 
surface (first impact at 306 sec past the atmospheric in- 
terface) bouncing more than 15 times before coming to 
a roll-stop. 

Pathfinder’s direct-entry strategy was in contrast to 
the Viking approach in which a propulsive orbit inser- 
tion was performed prior to descent and landing. As 
shown in Fig. 8, this decision significantly increased 
Pathfinder’s entry velocity, magnifying the design re- 
quirements on the entry, descent, and landing system. In 
addition, while the Viking landers were flown with an 
offset center-of-mass to achieve a lifting trajectory, Path- 
finder entered the Mars atmosphere at a near zero angle 
of attack to increase its drag force (reduce its ballistic 
coefficient) and simplify the design. In another simpli- 
fication over Viking’s active control strategy, 
Pathfinder’s flight at near zero angle of attack was main- 



Fig. 8: Mars Pathfinder and Viking atmospheric entry 
comparison. 

tained without propulsion. Instead, this angular orienta- 
tion was enforced through reliance on a 2 rpm spin rate 
and the inherent aerodynamic stability of the 70 degree 
sphere-cone aeroshell. 10 

ENTRY ANALYSIS 

Independent simulations of the Pathfinder entry, 
descent, and landing flight dynamics were developed at 
the Jet Propulsion Laboratory and the NASA Langley 
Research Center. At JPL, a three-degree-of-freedom 
translational equation of motion simulation was devel- 
oped with use of the Atmospheric Entry Program 
(AEP). 11 This simulation spanned the entire entry, de- 
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scent, and landing sequence (from the entry interface to 
landing). Additionally, the parachute descent and land- 
ing phases of the mission were analyzed with the 
Automated Dynamic Analysis of Mechanical Systems 
(ADAMS) multi-body simulation program in which both 
the translational and rotational equations of motion were 
solved. 12 At NASA-Langley, the six-degree-of-freedom 
version of the Program to Optimize Simulated Trajecto- 
ries 13 (6D POST) was used to analyze both the 
translational and rotational equations of motion from the 
atmospheric interface to parachute deployment. The 
three-degree-of- freedom version of POST was then used 
to simulate the parachute descent and landing phases of 
atmospheric flight. 

When coupled with the JPL orbit determination soft- 
ware (interplanetary cruise), this suite of atmospheric 
simulation tools provided a wide range of analysis capa- 
bilities. While the 6D POST and ADAMS simulations 
were believed to yield more accurate flight predictions 
(as a result of the more complete dynamic modeling), the 
three-degree-of-freedom AEP and POST simulations 
were used to provide an independent check of the higher 
fidelity results as well as rapid analysis capability. Good 
agreement in regard to the heati ng rate, parachute deploy- 
ment time, and other trajectory events was typical among 
this suite of Pathfinder-specific atmospheric analyses. 

The simulations relied on a consistent set of physi- 
cal models. 14 Atmospheric density and pressure profiles 
derived from Hubble Space Telescope and Earth-based 
microwave measurements of the Mars atmospheric tem- 
perature were employed. 1 ’ 5 Prior to entry, this model 
indicated that the atmosphere was likely to be signifi- 
cantly cooler with a lower dust content than the 
atmosphere observed by the Viking landers. Updates to 
this atmospheric model were performed during interplan- 
etary cruise as part of the entry operations procedure. A 
six-degree-of-freedom aerodynamic model developed 
from a combination of computational fluid dynamic cal- 
culations and existing wind-tunnel and ballistic-range 
data was employed. 10 This aerodynamic database was 
valid from entry (in the free molecular flow regime) to 
parachute deployment. After parachute deployment, 
aerodynamic drag predictions from Pioneer Aerospace 
(the parachute manufacturer) were used. 15 In the 
ADAMS multi-body simulation, additional aerodynamic 
relations were derived. 12 

The JPL and LaRC simulations were developed in 
the spacecraft design phase and utilized throughout the 
Pathfinder program in a Monte Carlo fashion to statisti- 
cally assess the impact of a range of off-nominal condi- 
tions on the flight system. 17 - 18 Outputs from these Monte- 


Carlo simulations were used in the design of the Path- 
finder heatshield, the entry, descent, and landing flight 
software, and to define numerous sets of hardware tests. 19 
In addition, these simulation results were imported into 
JPL’s spacecraft testbed and used to evaluate the per- 
formance of the entry, descent, and landing flight soft- 
ware for a range of off-nominal atmospheric flight con- 
ditions. 

In the operations mode, the 6D POST, 3D POST, 
and AEP simulations were used by the navigation team 
to update the estimation of the atmospheric flight condi- 
tions. Changes in Pathfinder's estimated atmospheric 
flight were expected as a result of improved estimates 
of the atmospheric model and atmospheric interface state 
vector as the spacecraft approached Mars. By modify- 
ing the entry, descent, and landing flight software pa- 
rameters, the operations team could inform the space- 
craft of its most likely atmospheric flight conditions. 
Without this update capability, the likelihood of a suc- 
cessful entry, decent, and landing (particularly, a suc- 
cessful parachute deployment) would have been ad- 
versely affected. As a result, these atmospheric flight 
simulations were employed at numerous opportunities 
preceding Pathfinder's entry, descent, and landing to 
update the onboard set of flight software parameters and 
estimate the landing site. This update process was one 
of the primary responsibilities of the operations naviga- 
tions team in the days prior to landing. 

ENTRY, DESCENT, AND LANDING FLIGHT 
SOFTWARE PARAMETERS 

Parachute Deployment Software Parameters 

Pathfinder’s entry, descent, and landing software was 
responsible for autonomously guiding the spacecraft 
from cruise-stage separation to a successful landing. 
Much of this sequence was timed by processes running 
in the flight computer, with events being triggered by 
pyrotechnic firings. Both a primary and backup system 
were developed to initiate deployment of the parachute. 
Without the luxury of an inertial measurement unit or 
gyroscopes, the primary parachute deployment algorithm 
relied solely on accelerometer readings and consultation 
of a pre-determined set of entry deceleration protiles, 
stored in curve-fit form. 20,21 This algorithm was initi- 
ated at cruise-stage separation and was responsible for 
interpreting the deceleration pulse, scheduling the ap- 
propriate time to initiate parachute deployment, and fir- 
ing the parachute mortar. 

The primary algorithm relied on initiation of a timer when 
the deceleration level reached 5 g. This is termed the g] 


5 



Flight Software: 

• Waits for deceleration pulse 
(accelerometers) 


Altitude, Deceleration • Computes time-to-go 

km Earth g’s • Schedules parachute deployment 



Fig. 9: Pathfinder flight software parachute deployment process. 


point. Twelve seconds after this timer initiation, a sec- 
ond deceleration measurement was taken (g 2 ). Based on 
the g 2 deceleration value, the onboard curve-fit was con- 
sulted to determine the time remaining until the para- 
chute could be safely deployed (t g0 ). This onboard pro- 
cess is illustrated in Fig. 9. As shown in Fig. 10, for the 
expected range of entry flight-path angles (-14.2 ± 1.0 
deg), a diverse set of deceleration profiles were expected, 
resulting in a relatively large range of parachute deploy- 
ment times. 



Time from gi point, sec 

Fig. 10: Range of deceleration profiles accommodated 
by parachute deployment algorithm. 

The objective of this in-flight software process was 
to deploy the parachute as close to a dynamic pressure 
of 600 N/m 2 as possible. Design requirements for the 
parachute included a dynamic pressure below 703 N/m 2 
and a Mach number greater than 1 .2. The 600 N/m 2 para- 
chute deployment target was determined iteratively 
through Monte-Carlo simulation, and maintained the 3- 


sigma high dynamic pressure below the 703 N/m 2 de- 
sign limit while providing sufficient time for the remain- 
ing elements of the descent and landing sequence to oc- 
cur nominally. Pre-launch Monte-Carlo simulations 18 
predicted a 99.1% probability of successful parachute 
deployment if the primary algorithm was successfully 
employed. 

The entry deceleration curve-fit used in this pro- 
cess was derived from a set of six-degree-of-freedom 
POST simulations specifically tuned to the latest entry 
state prediction and atmospheric model. The POST simu- 
lations were initiated with a dispersed set of entry states, 
terminating at 600 N/m 2 . These dispersed entry states 
were obtained by sampling the orbit determination co- 
variance at 0.1 deg entry flight-path angle increments. 
The maximum dispersion was defined by the three-sigma 
uncertainty in the orbit determination solution. As an 
example, the curve-fit loaded during spacecraft cruise is 
shown in Fig. 11. This linear curve-fit was estimated 
from a least-squares fit to the 6D POST data resulting in 
the following equation: 

t g0 = (40.14242 - g 2 )/ 0.28033 (1) 

This least-squares information was relayed to the 
spacecraft through two (g 2 4 g o) coordinates - a minimum 
(g2a) and maximum (g2b) expected value, 1 1 and 21 re- 
spectively. Using these values in the above equation, the 
minimum (t g0 b) and maximum (t goa ) may be computed 
as 68.284 and 103.956 sec respectively. In addition to 
defining the t g0 curve-fit, the g 2a and g 2 b values were 
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Fig. 11: Parachute deployment t g0 curve-fit onboard 
spacecraft prior to EDL parameter update #1 . 

also used to determine the validity of the accelerometer 
measurements during flight. As part of the fault protec- 
tion process, the primary algorithm was assumed to be 
invalid if the g 2 measurement was outside of the range 
of vales specified by g 2a and g 2b . Similarly, the com- 
puted t go was assumed to be invalid if it was outside the 
range specified t goa hy and t g0 b- 


ware parameters which determined the primary system’s 
curve-fit, fault-protection logic, and deceleration sam- 
pling strategy as well as the fixed-time backup could be 
updated during flight (prior to entry). Update criteria and 
command approval strategies were also established and 
simulated in several operations readiness tests. Param- 
eter updates were expected as the spacecraft’s entry state 
and predicted atmosphere varied. 

Predicted Landing Site Parameters 

The estimated landing site latitude and longitude 
were also loaded into the flight software as parameters 
which could be updated prior to encountering the Mars 
atmosphere. This landing site position estimate was used 
by the lander during surface operations to orient the high- 
gain antenna for Earth communications. The position 
accuracy required for successful high-gain antenna point- 
ing was approximately 1.0 deg (60 km). A parameter 
update to the landed estimate would have been required 
during cruise had the best-estimated landing site moved 
out of this 60 km tolerance. 

ENTRY OPERATIONS 


Due to the 8 hz flight software sampling frequency 
of the accelerometer measurements and potential data 
corruption, the gi and g 2 measurements may not occur 
precisely at 5 g and 12 seconds later. Hence, fault-toler- 
ant procedures were developed to estimate the true time 
at the gi point and the associated g 2 value 12 sec later. 
Windows, 1 .5 g and 1 .5 sec in duration, were established 
for the primary system to obtain valid data. If the sys- 
tem did not obtain a valid measurement in either of these 
windows, the backup parachute deployment algorithm 
would be enabled. Errors in the obtained measurements 
were minimized and the performance of the primary al- 
gorithm was enhanced by extrapolating from the true 
measurements to the desired times along pre-determined 
slopes (dtdg at the gi point and dgdt at the g 2 point). 


As shown in Fig. 12, the operations navigation func- 
tion was most critical in the 36 hours preceding entry, 
descent, and landing. In this period, four opportunities 
to update the entry, descent, and landing software pa- 
rameters and two opportunities to perform a contingency 
trajectory correction maneuver (TCM-5) existed. This 
final navigation function ended at approximately 6 am 
PDT on July 4 (entry -4 hours). 


July 2 I 
6 p.m. 12 a.m. 6 a.m. 


EDL 

Parameter 

Update 

Analysis 


♦ 

Update 

#1 ‘ 


TCM-5 

Analysis 


July 3 I July 4 

12 p.m. 6 p.m. 12 a.m. 6 a.m. 




♦; 

Update 

#2; 


!♦ ♦ 

Update Update 
j#3 #4 


E3 

♦ 

TCM-5a 


E3 

♦ 

TCM-5b 


12 p.m. 


EDL 

| 


If the flight software determined that accelerometer 
readings were invalid, the backup parachute deployment 
system would be enabled. This system initiates deploy- 
ment of the parachute at a fixed time, stored as a flight 
software parameter. 22 Pre-launch Monte-Carlo simula- 
tions 1 ^ predicted a 90% probability of successful para- 
chute deployment if the secondary system was employed 
and updates to this fixed-time strategy were performed 
as the spacecraft approached Mars. 

To minimize risk, parameter updates to the primary 
and secondary parachute deployment systems were de- 
signed into the operations navigation procedures. Soft- 


□ Orbit determination 
E3 Atmospheric flight analysis 
+ Command approval meeting 

Fig. 12: Operations navigation function timeline 
preceding atmospheric entry. 

EDL Update #1 

At approximately 8:00 pm PDT on July 2, 1 997 (en- 
try -38 hours) entry, descent, and landing update #1 was 
initiated. This ground analy si s process began by perform- 
ing orbit determination with the latest tracking data 
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included in the analysis. From this solution, a best-esti- 
mated trajectory was produced in which the entry 
flight-path angle was -1 3.942 deg. At this point in time, 
the entry covariance predicted a three-sigma flight-path 
angle uncertainty of ±0.8 deg. With this covariance a 
dispersed set of entry states were generated and six-de- 
gree-of-freedom atmospheric trajectory simulations were 
computed from the entry interface to the nominal para- 
chute deployment point (600 N/m 2 dynamic pressure). 
For each trajectory, the g 2 deceleration value and the t g0 
were computed, stored, and plotted. As shown in Fig. 13, 
a least-squares curve was then fit to this data for potential 
transmission to the spacecraft. At update #1, this linear 
least-squares fit produced the following equation; 

t g0 = (34.5 1997 - g 2 )l 0.2 1 802 (2) 



Fig . 13: Parachute deployment t g0 curve-fit loaded 
onto spacecraft as a result of EDL parameter 
update #7 process. 

The resulting values of t goa and t go b are 107,880 sec 
and 62.012 sec corresponding to a g 2 a of 1 1.0 g and a 
g 2 b of 2 1 .0 g, respectively. This data is plotted in Fig. 1 3 
which also presents the parachute deployment curve-fit 
which was loaded onboard the spacecraft at the time of 
this update opportunity. As shown in this figure, for the 
best-estimated flight-path angle at entry -38 hours, (a 
g 2 deceleration of 14.74 1 ), the onboard and current esti- 
mates of the nominal parachute deployment time differ 
by only 0.1 1 sec. In this case, such good agreement oc- 
curs because the best-estimated trajectory falls in the 
region where the two curve-fits fortuitously cross. For 
another entry angle, for example -14.06 deg where the 
g 2 deceleration value is 12.482, the two curve- fits dis- 
agree by approximately 2 sec with the update # 1 curve- fit 
implying a later parachute deployment (t g0 = 101.082). 

The best-estimated trajectory was then examined in 
detail such that values of dtdg and dtdg could be esti- 


mated ( 1 .660 and 0.745 respectively). Verification of the 
best-estimated 6D POST trajectory solution was per- 
formed through comparison with an AEP solution. The 
nominal parachute deployment time obtained from these 
two solutions differed by 1.95 sec (approximately 30 
N/m 2 in dynamic pressure) which was within the speci- 
fied tolerance. The backup parachute deployment time 
was then inferred from the 6D POST best-estimated tra- 
jectory (16:54:34.612 UTC, or 164.11 sec from the 
atmospheric interface). 

To determine if a parameter update was warranted, 
the current set of software parameters was compared to 
the onboard set as shown in Table 2. Pre-flight criteria 
had been established on the three most critical EDL pa- 
rameters. These general criteria stated that an EDL 
parameter update would occur if either the primary or 
backup parachute deployment time was predicted as being 
in error by more than 3 sec for the best-estimated trajec- 
tory or if the predicted landing site was off by more than 
1 .0 deg (60 km from the onboard value (for high-gain 
antenna pointing accuracy). As shown in Table 2, at entry 
-38 hours, without modification, the onboard backup para- 
chute deployment timer would have deployed almost 9 


Table 2. Entry, Descent , and Landing Software 
Parameter Update #/ 


Parachute 

Values 

Currently 

Difference 

Deployment 

Estimated in 

Loaded 

(Update 

Parameter 

Update Process 

Values 

Criterion) 

Primary algorithm 




parachute deployment 1 64. 1 2 

164.01 

0.11 (3.0) 

time, sec 




Fixed-time backup 
(UTC) 

16:54:34.612 

16:54:25.691 

8.92 (3.0) 

gi target 
deceleration, g 

5.0 

5.0 


gi to g 2 target 
interval, sec 

12.0 

12.0 


92a 

11.0 

11.0 


tgoa 

107.880 

103.956 


92b 

21.0 

21.0 


’gob 

62.012 

68.284 


dtdg @ g-j point 

1.66 

1.60 


dgdt @ g 2 point 

0.745 

0.84 


High-Gain 

Value 

Currently 


Antenna 

estimated 

loaded 


Parameter 

in update 

value 



process 



Latitude, 

18.82, 

19.2438, 

0.69° (1.0°) 

Longitude 

326.32 

326.9000 

41 .0 km 
(60 km) 
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sec early. Such a large error was a result of the more 
shallow entry angle, an update to the atmospheric model, 
and bit cutoff errors made in the onboard set of com- 
mands. As a result, a parameter update was deemed 
necessary. After a command approval meeting, the Table 
2 parameters were relayed to the Pathfinder spacecraft at 
on July 2, 1997 at approximately 1 1:00 pm PDT. 

While the ground analyses for EDL update pro- 
cesses 2-4 were all performed, the spacecraft remained 
close to its predicted path. As an example, as these op- 
portunities passed, the best-estimated entry flight-path 
angle was —13.942,-13.902, -13.896, and —13.914 deg 
at the four update opportunities. As a result, no further 
EDL software modifications were relayed to the space- 
craft. Hence, the EDL parameters presented in Table 2, 
Fig. 13, and equation 2 were the onboard set used dur- 
ing Pathfinder’s atmospheric flight. 

Landing Site Prediction and TCM-5 

As discussed previously, four trajectory correction 
maneuvers (TCM) were designed into the nominal Path- 
finder flight profile. As a result of planetary protection 
requirements, the first of these to actually place the space- 
craft on an intercept trajectory with Mars was TCM-3. 
After performing this 1 1.0 cm/s AV on May 6, 1997, 
orbit determination solutions indicated that the Pathfinder 
was on a trajectory with a steeper atmospheric entry flight 
path angle than desired. At this point, the predicted en- 
try flight-path angle was -14.84 deg; whereas, an entry 
flight-path angle of-14.2 deg was desired. As shown in 
Fig. 14, this resulted in a landing site prediction roughly 
equal in size, but offset to the northeast of the science 
requirement ellipse. Hence, prior to TCM-4, approxi- 
mately 60% of the predicted landing ellipse did not meet 



Fig . 14: Post TCM-3 landing site prediction. 


the mission requirements. Note that the ellipse shown in 
Fig. 14 represents a three-sigma prediction; that is, a 
99.7% probability existed that without further maneu- 
vers the final landing site would be within this predicted 
region. 

The fourth trajectory correction maneuver was ex- 
ecuted on June 25, 1997. This 1 .8 cm/s maneuver placed 
the spacecraft very close to its nominal path (best-esti- 
mated entry flight-path angle of -14.18 deg). 
Furthermore, after completion of this maneuver, the pre- 
dicted landing ellipse was completely contained within 
the science requirement ellipse measuring approximately 
180 x 70 km. This landing site prediction is shown in 
Fig. 15. 



Fig. 15: Post TCM-4 landing site prediction. 

As Pathfinder approached Mars, small changes in its 
predicted landing site were expected as a result of in- 
creased state knowledge (once the Mars gravity well w-as 
sensed) and Mars ephemeris errors. In the event that the 
predicted landing site drifted from the target ellipse, a 
contingency maneuver (TCM-5) would be decided upon 
at either entry -1 1 or entry -7.5 hours (see Fig. 12). 
Because this maneuver would occur so close to entry, a 
customized maneuver design could not be performed, (as 
was done for TCMs 1-4). Instead, in each of the TCM-5 
opportunities, the flight team had to decide whether to 
activate a command sequence from a pre-defined set al- 
ready loaded on the flight system. In this manner, through 
proper sequence selection, the predicted landing site 
could be moved in downtrack and crosstrack. 
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Figure 1 6 shows how the predicted landing site var- 
ied as the spacecraft approached Mars. Three-sigma 
ellipses at entry -24, entry -13, and entry -9 hours are 
shown. The last two ellipses represent the flight team’s 
knowledge at the time of the two TCM-5 decision win- 
dows. Each of the TCM-5 landed prediction ellipses is 
approximately 100 x 15 km. While the best-estimated 
trajectory was predicted to be more shallow (an entry 
flight-path angle of -13.9 deg) than just after TCM-4, 
the decision was made not to perform a TCM-5. This 
decision was a result of: (1) a majority of the TCM-5 
predicted landed ellipses lying within the science require- 
ment, (2) the science team not strongly objecting to a 
small overshoot of the desired ellipse, (3) the navigation 
team’s confidence that the predicted ellipses would sig- 
nificantly diminish in size (but not move in center) as 
the entry approached, (4) the small northeastern move- 
ment in the prediction ellipse from entry -1 3 to entry -9 
hours, and (5) the small, but finite risk associated with 
doing a propulsive maneuver so close to Mars arrival. 



Fig. 16: Landing site predictions in the 24 hours 
preceding atmospheric entry . 

As the spacecraft continued on toward Mars, the 
navigation team continued to refine its landing site pre- 
dictions. Just prior to entry, the size of the predicted 
three-sigma ellipse had diminished considerably as 
shown in Fig. 17. The best-estimated flight path angle 
was -14.06 deg. As a result, the navigation ellipse (ex- 
cluding atmospheric and aerodynamic effects) was 
centered on a landing site of 19.22 deg N latitude, 33.4 
degree W longitude. This ellipse was only 15 x 7 km in 



Fig. 1 7: Landing site predictions at atmospheric entry. 

size and was no longer oriented with its semi-major axis 
aligned to the flight direction. This axis orientation 
change is indicative of a small degree of uncertainty in 
the entry flight-path angle. 

The second ellipse shown in Fig. 17 represents the 
results of a Monte-Carlo analysis performed with the 
pre-entry orbit determination solution, including disper- 
sions in the atmospheric flight-path. This ellipse is 
approximately 40 x 15 km, with its semi-major axis bet- 
ter aligned along the flight direction. It is centered on a 
landing site of 19.15 deg N latitude, 33.51 degree W 
longitude. Prior to entry -24 hours, atmospheric disper- 
sions had been justifiably neglected in the navigation 
team’s landing site predictions. However, as shown in 
Fig. 17, as knowledge of the entry state increased, the 
relative significance of the atmospheric flight dispersions 
also increased. These dispersions, caused by aerody- 
namic and atmospheric uncertainty, more than doubled 
the size of the final landing ellipse prediction. 

The landing site determined by the Pathfinder sci- 
ence team once the spacecraft was safely on the Mars 
surface is also denoted in Fig. 17. This landed estimate 
was determined using the lander images to triangulate 
from observed surface features (craters, knobs, and 
peaks) as shown in Fig. 18. The science team landed 
estimate (19.33 deg N latitude, 33.55 deg W longitude) 
places the spacecraft just 27 km from the navigation tar- 
get, the center of the science requirement ellipse. A small 
discrepancy (approximately 5 km) between the final pre- 
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Fig. 18: Science team determination of Pathfinder landed site . 


entry landed predictions and post-flight position reck- 
oning exists. A majority of this error is believed to result 
from map-tie errors where the Mars surface features are 
not accurately reflected in an analytical latitude/longi- 
tude map. Additionally, a post-flight reconstruction 
analysis of the accelerometer data 5 has indicated that the 
spacecraft may have bounced/rolled as much as 1 km 
before the airbag system stopped. This theory has been 
corroborated by evidence produced by the Pathfinder 
imager which has located an object believed to be the 
Pathfinder backshell 1 .2 km southeast of the lander. Find- 
ing the backshell in this direction (an azimuth angle of 
143 deg) further reduces the difference between the pre- 
entry and post-flight landing site predictions. 

PARACHUTE DEPLOYMENT ALGORITHM 
SYSTEM PERFORMANCE 

Engineering reconstruction of the Pathfinder entry, 
descent, and landing system performance 5 is a critical legacy 
of this technology demonstration mission with significant 
implications for future exploration (e.g.. Mars 98 lander, 
Mars 01 lander, Mars 01 orbiter, etc.). Atmospheric flight 
data from two sets of accelerometers, the radar-altimeter, 
and pressure and temperature sensors were measured and 
have been returned from the spacecraft. 


Analysis of the accelerometer measurements con- 
clusively demonstrates that the parachute was deployed 
based on the primary algorithm. Post-flight processing 
of the recorded accelerometer data shows that the ve- 
hicle passed through 5 g at 58.29 sec after encountering 
the atmosphere. From the flight data, a g 2 deceleration 
of 1 2.482 g is observed, yielding a t g0 of 1 01 .08 sec and 
a parachute deployment 171.37 sec past the atmospheric 
interface. In contrast, the fixed-time backup parameter 
would have deployed the parachute at 164.1 1 sec (see 
discussion regarding Table 2). 

Figures 19 and 20 present the vehicle deceleration 
as a function of time from the atmospheric interface, taken 
directly from the accelerometer readings. These mea- 
surements are indicated in Fig. 20 from which the 32 Hz 
sampling rate is evident. The downward spike in the entry 
vehicle deceleration, seen between the accelerometer 
samples at times 171.375 and 171.406, marks the firing 
of the parachute mortar. This temporary but significant 
acceleration (decrease in deceleration) is a result of the 
downward force imparted to the entry vehicle as the para- 
chute mortar fired upward through the backshell. Para- 
chute inflation required roughly 1.25 sec, with a maxi- 
mum parachute snatch load of approximately 6.5 g. Pre- 
liminary indications from the trajectory reconstruction 
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SUMMARY 

On July 4, 1997, the Pathfinder spacecraft success- 
fully landed on the surface of Mars. After a complex 
entry, descent, and landing process, the system landed 
just 27 km from its target point. In the present paper, the 
atmospheric entry and approach navigation activities re- 
quired in support of this mission have been discussed. 
In particular, the flight software parameter update and 
landing site prediction processes performed by the Path- 
finder operations navigation team were described. 


Fig. 19: Parachute deployment dynamics extracted 
from Pathfinder flight data. 



Time from atmospheric interlace, sec 

Fig. 20: Time of parachute mortar firing extracted 
from Pathfinder flight data. 

analysis indicate that the parachute deployment altitude 
was approximately 7.9 km above the Mars reference el- 
lipsoid (9.5 km above the surface) at a dynamic pressure 
of approximately 588 N/m 2 and a Mach number of 1 .7 1 . 
This is extremely close to the parachute deployment de- 
sign target; hence, the primary parachute deployment 
algorithm is believed to have performed well. 

It is also interesting to note that the original fixed- 
time backup parachute deployment time, prior to param- 
eter update process # 1 , would have resulted in a deployed 
parachute 16.2 sec early at a dynamic pressure close to 
850 N/m 2 (well above the design limit of 703 N/m 2 ). 
However, after modification during parameter update #1 , 
the backup fixed-time deployment would have occurred 
7.3 sec early at a dynamic pressure of approximately 
695 N/m 2 . This demonstrates that the parameter update 
process performed by the operations navigation team in 
the final days of interplanetary cruise increased the prob- 
ability of mission success. 


The final set of flight software parameters was relayed to 
the Pathfinder spacecraft following the first parameter 
update process on July 2, at 1 1 :00 pm PDT. This change 
set included modification of parameters affecting both 
the primary andbackup parachute deployment algorithms 
Although EDL parameter update 2-4 were all performed, 
no further updates were deemed necessary by the opera- 
tions team. Preliminary evaluation of the accelerometer 
flight data indicates that the Pathfinder parachute deploy 
mortar was activated by the primary software algorithm 
at the appropriate conditions (588 N/m 2 and a Mach 
number of 1.71). This performance validates the design 
of the parachute deployment software algorithm. F urther- 
more, post-flight analysis has demonstrated that without 
the software parameter update and on the backup timer, 
the parachute deployment design constraints would have 
been significantly exceeded (a dynamic pressure close to 
850 N/m 2 ). Hence, the probability of mission success was 
increased through the update process. 

As the spacecraft approached Mars, the operations navi- 
gation team continued to refine its landing site predic- 
tions. At the two TCM-5 decision opportunities the pre- 
dicted landing ellipse was approximately 100 x 15 km 
in size and largely within the science requirement re- 
gion. As a result, this contingency maneuver was not 
performed. As the spacecraft continued its Mars ap- 
proach, the size and orientation of the predicted landed 
ellipse dramatically changed. At the time of entry, the 
predicted three-sigma landed ellipse size was 15x7 km, 
neglecting aerodynamic and atmospheric uncertainty. 

At entry, the three-sigma flight-path angle uncertainty 
was 10.05 deg, centered on a nominal value of -14.06 
deg. With such a small uncertainty on the entry state, 
the effect of aerodynamic and atmospheric uncertain- 
ties on the flight path was shown to be significant. With 
these uncertainties included, a six-degree-of-freedom 
Monte-Carlo analysis resulted in a three-sigma landed 
ellipse estimate of 40 x 15 km. The center of this esti- 
mate is within 0.5 km in downtrack position of the post- 
flight landing estimate produced by the Pathfinder sci- 
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ence team (with a crosstrack difference of approximately 
5 km), demonstrating the importance of atmospheric 
flight modeling to precision-landing site prediction. The 
remaining crosstrack difference is postulated to be a re- 
sult of map-tie errors in which the Mars surface features 
are not accurately represented in the analytic map. 
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